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We study local density of electron states of a two-dimentional conductor with a smooth disorder 
potential in a non-quantizing magnetic field, which does not cause the standart de Haas-van Alphen 
oscillations, ft is found, that despite the influence of such "classical" magnetic field on the average 
electron density of states (DOS) is negligibly small, it does produce a significant effect on the DOS 
correlations. The corresponding correlation function exhibits oscillations with the characteristic 
period of cyclotron quantum huic- 
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I. INTRODUCTION 

In the clean homogeneous electron gas the wave func- 
tions of electrons are plane waves and the density of elec- 
tron gas is constant in space. In disordered conductors 
electrons are scattered by impurities, which change their 
wave functions from the plane waves. This, in turn, re- 
sults in spatial variations of the electron density. It is 
appropriate to describe those variations by introducing 
the local density of electron states, i^(e,r), which is de- 
termined by the equation 



a 



(1) 



where index a specifies the electronic states, and the fac- 
tor of two reflects the spin degeneracy. The distribution 
function of the local DOS at the fixed energy in open 
metallic disordered samples was studied in many papers 
(see e.g. with the emphasis on the rare, nontypical 
fluctuations. It was found that although the local DOS 
distribution is close to the Gaussian one, it has slowly de- 
caying logarithmically normal asymptotics. Prigodin |^] 
studied correlation function of the density of the electron 
states of a two-dimensional system at different energies 
in relation to the NMR line shape. 

It is well-known that strong magnetic field modifies 
the single-particle density of electron states, both local 
and average, due to the Landau quantization. In a two- 
dimensional electron gas the quantization leads to a peak 
structure in the average density of states, which is re- 
vealed in tunneling experiments as peaks in the depen- 
dence of the tunneling conductance on the applied bias, 
see, e.g., Ref. The form and width of these peaks are 
determined Q by the disorder. 

In a weak magnetic field the distance between the Lan- 
dau levels, hujc, is smaller than their disorder-induced 
width. As a result, in such "classical" magnetic field, 
oscillations in the average density of states caused by 
the Landau quantization become exponentially small Q , 
cx exp[— 27r/(ti;c'^s)]- Here Ts is a quantum lifetime of an 
electron. 



The goal of the present paper is to show that, despite 
such "classical" magnetic field does not influence the av- 
erage DOS, it does produce a significant effect on the 
correlation function of the local density of states fluctua- 
tions 
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Here 5i'{e,r) — i^(e,r) — i^q is the local deviation of the 
DOS in point r from its average value, i^q — m/Trh^ , m is 
the electron mass, and brackets ( ) denote averaging over 
the random impurity potential. Clearly, the correlation 
function depends only on difference of energy agruments, 
P(ei,e2,r) = P(ei-e2,r). 

The effect of the classical magnetic field on the DOS 
correlation function becomes pronounced if the disorder 
potential has a correlation length much larger than the 
Fermi wave length. In such a potential, electrons expe- 
rience small-angle scattering, and their transport relax- 
ation time Ttr, is much larger than r^. Thus there exists 
a range of magnetic fields, in which Landau quantization 
is suppressed {uJcTs ^ 1), while classical electron tra- 
jectories are strongly affected by the field {uicTtr ^ !)■ 
In this regime the correlation function, P(ei — 62) is 
strongly enhanced with respect to the zero magnetic 
field case and exhibits peaks as a function of energy dif- 
ference ei — £2 with the distance between peaks equal 
to the cyclotron quantum, HcOc. For the macroscopi- 
cally homogeneous sample the shape of the n-th peak, 
|(ei — £2) — Tihu!c\ ^ ?iWc/2, in the local DOS correlation 
function is given by: 



where 
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and Ep is the Fermi energy. As n gets bigger, the width 
of the peaks increases and their height decreases, so that 
eventually the oscillatory structure is washed out. The 
total number of resolved peaks is of the order of ^Wc'^'tr- 
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Sensitivity of the correlation function P to the classical 
magnetic field comes from the fact, that this function is 
directly associated with the self-crossing of classical elec- 
tron trajectories. We denote the probability for an elec- 
tron to complete a loop of self-crossing trajectory over 
time t as K{t). The correlation function, P(ei — £2), 
turns out to be proportional to the Fourier transform of 
this return probability, 

/■OO 

Jo 

The strong enough, cjc^tr ^ 1, magnetic field curves the 
electron trajectories, significantly affects the return prob- 
ability and, in turn, leads to specific correlations in the 
local DOS at ei — £2 ~ ^^c- 

For long time scales t ^ rtr, the function K{t) can 
be found from the diffusion equation. It gives K{t) oc 
{Dt)~^ for the two-dimensional case {D is the diffusion 
coefficient). The Fourier transform, K{u>), is propor- 
tional to ln(w), which leads to the well-known log- 
arithmic form of the local DOS correlation function with 
the renormalized by the magnetic field diffusion coeffi- 
cient. 

At short time scales, i ^ rtr, electrons move ballisti- 
cally along the cyclotron orbits. Provided that oJcTtv ^ 1, 
during the time t electron may return to the initial point 
many times. Multiple periodic returns of electron pro- 
duce peaks in the probability Fourier transform K{uj) at 
energies, which are multiples of the cyclotron quantum. 
Correlation function P{ei — £2) oscillates with the same 
period, which is reflected by Eq. (^). 



P(£i -e2,r) - j—^ [2Re(t;iJ(r,r)C?,^(r,r)) 

+ <e«(r,r)e«(r,r))+(e,^(r,r)e,^(r,r))] . 

(7) 

The averages of the type {G^G^} and {G^G^} can be 
neglected as they do not contain contributions associ- 
ated with the electron trajectories longer than Xp and, 
thus, do not produce energy dependence of P(ei — £2, r) at 
£1 — £2 <C Ep. To the contrary, averag ed product (G^G"^) 
is determined by long electron trajectories (see e.g. Ref. 
Q). In general, the product of two exact Green's func- 
tions t/^ (ri, Y2)Gf'^ (fs, Y4) oscillates rapidly with the dis- 
tance between its arguments, so that the function 

iC(£i,£2,ri,r2,r3,r4) = (^i^ (ri, r2)e,^, (rg, r4)) (8) 

averages out. This is no longer the case if its argu- 
ments are close to each other pairwise. Namely, the sizes 
|ri - r4|, |r2 - ral or, alternatively, |ri -v^\, |r2 - r4| of 
spatial domains defining the ends of a trajectory should 
be small enough (less than vpTs) so that electron propa- 
gation in these two domains could be described by plane 
waves. 

If the ends of trajectories are separated by the distance 
exceeding electron wave-length, |ri — r2| Af, one can 
relate the function K to the generalized classical correla- 
tion functions - diffuson KP and Cooperon KP , which are 
given by the sum of all ladder and all maximally-crossed 
diagrams respectively. Namely, 



II. DERIVATION OF THE DOS CORRELATION 
FUNCTION 

Now we derive expression for the correlation function 
of the local DOS, P, valid for arbitrary electron energies. 
We omit the Planck constant in all the intermediate for- 
mulas. The DOS, Eq. (|), can be rewritten in terms of 
the exact retarded and advanced Green's functions of an 
electron in the following way: 
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K£,r) = -[c;,^(r,r)-C;f(r,r) 



where 



C?f(ri,r2)=5] 



C(l-l)V'a(r2) 
£ — £q -f iO 



(5) 
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and C/^(ri,r2) = [C/f^(r2, ri)]*. Single electron wave 
functions ipa (r) satisfy the Schrodinger equation for non- 
interacting electrons, Hgipa = (fia + Ep)^a, where Hq = 
— {h^ /2m)V'^ + Ur(r), and Ur{r) is the random potential. 

With the help of Eq. (||), the DOS correlation func- 
tion, Eq. (^), can be rewritten in terms of the ensemble- 
averaged products of the electron Green's functions: 



X(£i,£2,ri,r2,r3,r4) 

|ri - r4|, |r2 - ra] < vpTs 



K(£i,£2,ri,r2,r3,r4) 
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Here ~ ppUi, where = (cos 0^, sin^j) is a unit ver- 
tor with the direction determined by the angle 4>i. 

In the opposite limit, when the all four points ri, r2, 
r3, and r4 coinside, both ladder and maximally-crossed 
diagrams contribute to Eq. (^. As a result, the DOS 
correlation function, Eq. (^, contains both the diffuson 
and the Cooperon contributions: 
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P(ei,e2,r) = Re (r, r) + Ce,-., (r, r)] . (11) 

TTVq 

Here V and C are the diffuson, KP and the Cooperon, 
averaged over the initial and the final directions of 
the electron momentum: 



ri,r2 



27r 27r <^i-'^ 



(ri,0i;r2,</.2), (12) 
^(ri,0i;r2,(/)2). (13) 



As one sees, calculation of the DOS correlation func- 
tion reduces to the analysis of two classical correlation 
functions, T) and C. Provided that we are interested in 
the DOS correlation function in the presence of the mag- 
netic field, the problem can be further simplified. Indeed, 
as it is well-known, the diffusion and the Cooperon de- 
pend quite differently on the magnetic field (see e. g. 
Ref. 1^). In particular, C is exponentially suppressed if 
the magnetic length, \h = \J cTi/eH^ becomes smaller 
than the transport relaxation length Itr- We, in fact, 
assumed a much stronger condition, lJcTit S> 1, for the 
magnetic field. Thus the Cooperon term in Eq. ( pi] ) can 
be neglected in our case. On the other hand, the diffu- 
sion term in Eq. (|l^) is meaningful and will be analyzed 
below. 



III. DOS CORRELATION FUNCTION FOR AN 
INFINITE TWO-DIMENSIONAL ELECTRON 
GAS 

Let us first calculate the local DOS correlation func- 
tion, P(ei — £2), in the macroscopically homogeneous 
sample. The generalized diffuson, /Cj'(ri, (j)i;r2, (f>2), sat- 
isfies the Boltzmann equation (see e.g. Ref. ^) describ- 
ing the scattering of electrons on impurities in the pres- 
ence of the magnetic field. In the special case Tg ^ Ttr 
we are interested in, small angles scattering dominates 
the collision integral. With account for this simplifica- 
tion, the transport equation for /Cj^(ri, i/ii; r2, (^2) takes 
the Fokker-Planck form: 
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Equation (|Tj) describes electron motion along the cy- 
clotron orbit accompanied by the angular diffusion 
caused by scattering on a random potential. The solution 
of this equation will give us the Fourier transform. 



/•oo 

C(i;2)=/ /C^(t,i;2)e 
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of probability density, /C'''(i, ri, r2, ^2), for electron 
which starts at moment t = in point ri with the direc- 
tion of momentum (pi to arrive at moment t to the point 
V2 with momentum direction (j)2- 



In order to solve Eq. (|lj) it is convenient to introduce 
new spatial variables which correspond to the center of 
the electron cyclotron orbit: 



R = r + Rc[n x z]. 



(15) 



Here i?c is a cyclotron radius, and z is a unit vector paral- 
lel to the magnetic field. Changing variables in Eq. (|lj), 
and performing the Fourier transformation from R2 to 
q, we obtain: 
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2 R, \"^'^d^2 502 
X /C^(Ri,(/.i;q,</.2) = 27r5(0i-02)e-^'i^^ (16) 

We seek for the solution of Eq. (^6|) in the following 
form: 



/C^ (Ri , ; q, 02 ) = ^ e"*^ F„ (.., q; 



(17) 



Substitution of Eq. ( pT| ) to Eq. ( |T^ ) results in a linear sys- 
tem of equations for F„. At small enough wave vectors, 
qRc <C w^Ttr/dwl -I- ujc), terms in the square brackets in 
the l.h.s. of Eq. (|6|) become small, the equations cor- 
responding to different n become independent, and we 
obtain a solution for Fn{LO, q; 0i) in the form: 



K(w,q;0i) 



g— in0ig-iqRi 



-i(ti) — nUJc) + 
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The inverse transformation of variables immediately 
yields now the solution of Eq. ([T^: 



/CS'(ri,0i;r2,02)= / J^^^'""^''-'' 



gjn(02-0l)gj-Req[(n2-ni)xz] 



-i{uj — nuJc) + 
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After substitution of Eq. ( pj| ) into the Eqs. (11) and 
(|l^), and subsequent integration over angles, we obtain 
the following expression for the correlation function of the 
local DOS of a homogeneous two-dimensional conductor 
in the classical magnetic field: 



P{ei -62) = 



rfq 



(2^: 



-RePei-e2(q)> 
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(20) 
(21) 
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Here J„(z) is a Bessel function. At small frequencies, 
uj = ei — €2 't^ h/rtr, the n = term in Eq. domi- 
nates, and 



-iu; + R^^qy2Tt,' 



Overall energy dependence of the correlation function 
of the local density of states for an infinite sample, ob- 
tained by numerical analysis of Eqs. (pO|)-(pl|) is pre- 
sented in Fig. ^ The DOS correlation function exhibits 
strong oscillations with the period close to hujc- 
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FIG. 1. Energy dependence of the local density of states 
correlation function, P(ei —£2), of the macroscopically homo- 
geneous sample in the classical magnetic field, obtained by 
numerical analysis of Eqs. (|2o|)-(^l|). Parameter ujcTtr = 10. 

This limit corresponds to the diffusion regime with the 
diffusion coefRcient D ~ R^/2Ttr rcnormalized by the 
magnetic field. The correlation function P(ei — 62) de- 
pends logarithmically on ei — £2 in this limit: 



P{£l - £2) TT-F^ In 



ei - £2 Ttr 



(22) 



At large frequencies, w = ei — e2 S> h/rtr the correla- 
tion function P(ei — £2) exhibits peaks at ei — 62 close to 
multiples of cyclotron quantum, nhujc- The form of the 
n-th peak is given by: 



P(ei - e2) 



2.^2 



nEpTtr 



Re {In[a{n, Suj)]Kn[a{n, 5uj)]} , 



(23) 



where hSto = ei — e2 — nhujc was introduced instead of 
ei — £2, and a{n,Suj) = 2(n^ — iStuTtr). Functions In{a) 
and Kn{a) are the modified Bessel functions. For large 
a we can use the asymptotical relation, In{a)Kn{a) ~ 
l/2a, and arrive at the resulting Eq. (^) that describes 
energy dependence of the local DOS correlation function 
in the vicinity of n— th peak. 



IV. OSCILLATIONS OF THE DOS FOR 
TUNNELING INTO THE EDGE OF A 
TWO-DIMENSIONAL ELECTRON GAS 



The tunneling density of states i^(e, r) is directly re- 
lated to the tunneling differential conductance G{V) of 
a point contact attached to the two-dimensional gas, 
iy{eV,r)/L'o = G{V)/Go (here Go is the average linear 
conductance at zero magnetic field) . Thus measuring the 
conductance correlation function {SG{V)SG{V + AF)) 
one can determine the DOS correlation function, 
P{eAV) = (SG{V)SG{V + AV)) /Gl Here SG{V) = 
G{V) — Gq. The tunneling DOS we studied so far is re- 
lated to tunneling into the "bulk" of a two-dimensional 
electron gas, see Fig. For GaAs heterostructures, 
however, there exists a well developed method of form- 
ing point contacts for lateral tunneling into the edge of 
a two-dimensional electron gas (see, e.g., the review of 
Beenakker and van Houten, Ref. 0). The edge affects 
electron trajectories and thus alters the correlation func- 
tion of the tunneling density of states. Below we esti- 
mate P{ei — £2) for the specific case of lateral tunneling, 
schematically shown in Fig. ^b. We demonstrate that 
the oscillatory pattern of P at energies larger than fiLUc 
persists, although the amplitude of oscillations becomes 
smaller than in the case of tunneling into the bulk. 



tunnel contact 




b) 




tunnel contact 



FIG. 2. Two possible tunneling experiment that enable 
to measure properties of the tunneling density of electron 
states: (a) the point-like tunnel contact is attached to a 
two-dimensional conductor far from its edges; (b) tunneling 
occurs at the edge of the two-dimensional electron gas. 
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In order to find the conductance correlation function, 
one should, according to Eq. (|l|) , find the Fourier trans- 
form of the return probability, V^t), for an electron emit- 
ted from the contact right at the edge of the electron gas. 



case = 0, which correspond to the center of electron 
cyclotron orbit having the initial coordinates Rx — 0, 
Ry — Re- Our goal now is to obtain probability to find 
the center of orbit again in the same point after time t. 
During time t center of orbit diffuses in vertical (see Fig. 
S) direction on a distance 



Ar = 2R sind) 

X c 




Tunnel contact 



FIG. 3. Drift of an electron in the magnetic field caused 
by the multiple specular reflections from the boundary of the 
two-dimensional electron gas. 

Let us consider an electron which is emitted from a 
point-like tunnel contact attached to the edge of the 
two-dimensional conductor at a moment t = with the 
initial velocity characterized by the angle For any 
nonzero 4>i electron experiences multiple reflections from 
the boundary of the two-dimensional electron gas (see 
Fig. H). The boundary of the electron gas is usually 
smooth, so that we assume this scattering to be purely 
specular. Those multiple scattering events lead to a drift 
of the guiding center of electron orbit along the boundary 
of the 2DEG with the velocity Vd, which is given by: 



V2 



2tt 



-Vp 



Rr. Rii 



Rc 



for Ry less than Rc and is zero otherwise. Here R^ 



(24) 



and R,, 



are the coordinate of the center of electron orbit, which 
in the initial moment t = are: 



Rx = Rc sin 01 , Ry — Rc cos , 



(25) 



and ARx is defined in Fig. |3| In the absence of disor- 
der, drift prevents electron from return to the contact, 
and the return probability, I?(t) = at t > 0. Disor- 
der, however, makes the return probability nonzero. In 
fact, disorder leads to two effects: (1) motion along the 
cyclotron orbit is accompanied by the angular diffusion, 
and (2) in addition to the boundary-induced drift, the 
guiding center of the electron cyclotron orbit diffuses in 
the direction perpendicular to the boundary. As we will 
see, for small enough initial angles (f>i these two effects 
can, in fact, overcome the boundary-induced drift of elec- 
tron away from the contact. 

From Eqs. (|2^) and (^5|) we see that the larger ini- 
tial angle 01 is, the faster electron drifts away from the 
contact. In the view of this fact, let us start from the 



AY{t) = \Ry{t)-Rc 



'Dt, 



(26) 



where D = Rc/2Ttr is a diffusion coefficient. During the 
same time interval t, the center of orbit will travel along 
horisontal axis on a distance 



AX 



t pt 

Vd{t')dt' - VF 
^0 




AY{t'] 



Rc 



dt' 



^5/4 



.1/4 ■ 



(27) 



Here we exploited Eqs. (|24D and (p6[). One sees that the 
probability to find the center of electron orbit in the ini- 
tial point after time t decreases rapidly with time, 



K{t) 



1 



3/4 



cx 



AX Ay RcVptV^ ■ 



As the result, in the presence of a boundary, contribu- 
tions to the electron return probability coming from the 
trajectories with two and more revolutions along the cy- 
clotron orbit are small and can be neglected, while the 
main contribution comes from trajectories which involve 
only one revolution between start at t — and finish at 
the moment t ^ tc = Itx jiOc- Let us study now this latter 
contribution. 

If there were no boundary, the probability density 'D{t) 
for the electron to return to the initial point at time 
t = tc^ bt, where St tc, could be easily obtained from 
the solution, Eq. ([l9|), of the transport equation ([l^). 
In fact, one puts ri = r2 in Eq. (p^), and integrate it 
over all possible values of 0i and 02 taking into account 
that we are interested in the trajectories which are close 
to a single cyclotron loop. As the result we obtain the 
return probability density which has a strong maximum 
a.t t — tc with the amplitude depending on the amount 
of disorder in the system: 



V{t « tc) 



V2t 



'iOcTtr 



Rl 



■ exp 



t-tc 
tc 



(28) 



This equation is valid in the absence of the boundary, i.e. 
for a homogeneous system. Clearly, for such system tra- 
jectories with different initial angles 0i contribute equally 
to the Eq. (28). For the system with the boundary this 
is obviously not the case. Namely, only a small fraction 
of trajectrories with 0i ^ 1/ ^/uJcTtr contribute, for which 
the disorder-induced uncertainty of electron position ex- 
ceeds the shift ARx = Vdtc, see Fig. |. Thus in the 
presence of the boundary the return probability density. 
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'Dh{t), can be estimated by multiplying Eq. ( p8| ) by a 
small factor 1 / ^lOcTtr'- 



■Db{t) 



(29) 



According to Eqs. (|l l|) , the correlation function of the 
DOS at the edge of the two-dimensional electron gas, 
P(ei — £2), is determined by the Fourier transform of 
I?f,(i) given by Eqs. (p8[)-(p9|). Performing the Fourier 
transformation, we finally obtain: 



P(ei-e2) 



■ cos ZTT 



X exp 



27r 

LUcTtr 



flLOr 



ei - £2 



TlLJr 



(30) 



One sees that the correlation function exhibits harmonic 
oscillations with the period fiujc up to the energies of the 
order of ei — e2 ~ huj ^^uJ^tu- . The amplitude of these 
oscillation is iOcTtr 3> 1 times smaller than in the case of 
vertical tunneling into the bulk of the two-dimensional 
electron gas, see Eq. (0). 



V. DOS CORRELATION FUNCTION IN AN 
INTERACTING SYSTEM 



In order to observe experimentally the oscillations of 
the DOS correlation function predicted in the present pa- 
per and given by Eqs.(||) and (^0|), one has to distinguish 
them from the mteraciion- induced oscillations of the den- 
sity of the electron states. The easiest way to do this is 
to fix two of the arguments of the correlation function. 



t\ — e^' and ti — , and then measure P as a function 
of the shift in the chemical potential, e^'' — . 



,(2 



VI. CONCLUSIONS 

In summary, we study properties of the two- 
dimentional conductor with a smooth disorder potential 
in a magnetic field. It is known that the average den- 
sity of states of such a conductor is hardly modified by 
the magnetic field \&vlv^ (x exp(— 27r/wcTs)] as long as 
L^cTs ^ 1- We show that despite such "classical" mag- 
netic field does not influence the average DOS of the con- 
ductor, it does affect strongly the correlation function 
of the local density of states, P(ei — £2). Namely, pro- 
vided that ijJcTtr ^ 1, the correlation function P(ei — £2) 
aquires an oscillatory structure with the characteristic 
period hiOc- This structure can be observed in tunnel- 
ing experiments on both vertical tunneling into the bulk 
of the two-dimensional conductor, and lateral tunneling 
into the edge of the conductor. 



Until now we have completely disregarded effects of 
the electron-electron interaction. It is known, however, 
that this interaction has a crucial effect |^ on the tun- 
neling DOS of the disordered conductor. Namely, inter- 
action leads to a strong energy dependence of the single- 
particle density of electron states for the energies close to 
the Fermi level. As a result, the density of states must 
be written as a function depending both on the position 
of the Fermi level and on the electron energy measured 
]rom the Fermi level: 
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(31) 



In the two-dimensional system v{e — ep^ep) has a loga- 
rithmical singularity |p| at small e — ep <^ Ti/rtr^ and can 
be quite pronounced p.liol even at large e — S> 'h/rtr- 
In particular, in the classical magnetic field v{t — ep, ep) 
is an oscillating function [p^ of e — with a character- 
istic period of cyclotron quantum hujc. 

As a consequence of Eq. (|3l|), the DOS correlation 
function for an interacting system is a function of three 
arguments: 



P(ei,e2) = P(ei-e('^ - -^'^ 



"P 7 ^2 ^p 



J2) (2)- 



(32) 
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